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The theory of quantum trajectories is applied to simulate the effects of quantum noise sources 
induced by the environment on quantum information protocols. We study two models that generalize 
single qubit noise channels like amplitude damping and phase flip to the many-qubit situation. 
We calculate the fldelity of quantum information transmission through a chaotic channel using the 
teleportation scheme with different environments. In this example, we analyze the role played by the 
kind of collective noise suffered by the quantum processor during its operation. We also investigate 
the stability of a quantum algorithm simulating the quantum dynamics of a paradigmatic model of 
chaos, the baker's map. Our results demonstrate that, using the quantum trajectories approach, we 
are able to simulate quantum protocols in the presence of noise and with large system sizes of more 
than 20 qubits. 
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I. INTRODUCTION 

It would be highly desirable to implement quantum 
protocols using processors perfectly isolated from the en- 
vironment, since this is one of the main sources of error 
in quantum computation (there are also system specific 
imperfections, but here we will only address the envi- 
ronmental problem). Unfortunately, this is not possible. 
Quantum hardware will naturally become entangled with 
the environment during its operation. Thus, if any hope 
of profiting from the benefits of quantum computation is 
to be kept, understanding and controlling quantum noise 
effects is essential. On the other hand, the study of open 
systems is of interest in several fields, from both theoret- 
ical and experimental points of view P, 0, 0, ■ 

Factoring large integers in polynomial time has been 
the milestone discovery that set out quantum computa- 
tion as a major research topic Nevertheless, in the 
short term, few qubits quantum computers, this kind of 
calculations will be necessarily out of reach, since they 
involve large systems. Then, it is reasonable to focus 
on understanding the behavior of accessible first realiza- 
tions. It is interesting to remark that, with a few tens of 
qubits, quantum simulations of systems studied in quan- 
tum chaos (like the quantum baker's map &\, the quan- 
tum kicked rotator Q, and the quantum sawtooth map 
@) would outperform any calculation that can be done 
with present day supercomputers. A first step in the 
simulation of quantum chaos models has been the imple- 
mentation of the quantum baker's map on a three-qubit 
NMR-based quantum processor 0. 

With this situation in mind, it is natural to ask to what 
extent we can know and control the operability and the 
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stability of a quantum computer of this size. Theoretical 
studies rely on evaluating the reduced density matrix of 
the system (obtained after tracing out the environment), 
often in terms of various different approximations. It 
would therefore be desirable to give an answer for generic 
quantum protocols and noise models, using exact calcu- 
lations. In this work, we propose the numerical simu- 
lation of superoperators as a way to do it. By means 
of quantum trajectories techniques we can reach results 
for system sizes for which the implementation of chaotic 
maps becomes relevant for theoretical studies in the field 
of quantum chaos, having the chance to include any kind 
of environmental effects. 

Instead of solving the density matrix directly, quan- 
tum trajectories stochastically evolve the state vector of 
the system, and after averaging over many runs the same 
results for the outcomes of any observable are obtained. 
The use of quantum trajectories in the field of quantum 
information has been pioneered by [lol [Tll | . In Ref. ^3 ■, 
we applied this formalism to study the effect of a dissipa- 
tive environment on the quantum teleportation protocol 
[Tsj through a large chain of qubits. This situation mod- 
els the transmission of quantum information through a 
chaotic quantum channel. Here we review this model and 
present new results for different noise channels. 

As mentioned before, one of the main near future ap- 
plications of quantum computation is the simulation of 
dynamical systems that are of great interest in quantum 
chaos. In this work we focus on the quantum baker's 
map, one of the most important examples of this kind 
of systems. This map is fully chaotic and its quantized 
version consists of conveniently selected quantum Fourier 
transforms. We model the environment through a phase 
fiip channel and study the fidelity of quantum computing 
of the quantum baker's map in this noisy environment. 
We note that the fidelity has already been computed ex- 
perimentally with a three-qubit NMR quantum processor 
[9j. Hence, for the design and construction of quantum 
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hardware with a larger number of qubits, simulations like 
those performed in this paper will become essential. 

This paper is organized as follows. In Sec. ^ we 
present a brief explanation of the theory of quantum tra- 
jectories and connect it with the quantum operations ap- 
proach to the density matrix evolution. We also refer to 
the master equation formulation. Afterward, in Sec. 11111 
we make a short description of the noise channels that 
we use in the calculations, paying special attention to 
the amplitude damping models. In Sec. IIVI we review 
the quantum teleportation protocol presented in |l2l | , fo- 
cusing on the different dissipative processes. In Sec. 
we study the behavior of the fidelity of the quantum al- 
gorithm for the baker's map in the presence of a phase 
flip noise. Finally, in Sec. lVIl we present our conclusions 
and outlook. 



II. MASTER EQUATION, SUPEROPERATORS, 
AND QUANTUM TRAJECTORIES 

There is a close relationship among the master equa- 
tion in Lindblad form, the superoperator formalism and 
the quantum trajectories theory. It is useful to review 
this relationship and to motivate the use of the quan- 
tum trajectory approach when studying open systems 
like quantum processors. We start with the special case 
of the master equation formulation of the problem and re- 
late it with simulations using quantum trajectories tech- 
niques. Then, in Sec. Ill Bl we generalize this establishing 
the connection with the broader quantum operators for- 
malism. 



A. Master equation and quantum trajectories 

Real systems interact with the environment, and, as 
mentioned before, these interactions are usually referred 
to as quantum noise |14L Il5| . Several models account 
for different sorts of system-environment interactions, the 
particular choice depending on the nature of the system 
and environment under consideration. Such open quan- 
tum systems in general cannot be described by a pure 
state, but rather by a mixed state. Their evolution takes 
density matrices to density matrices. This allows the 
evolution from system's pure states to mixed ones, and 
also, in some noise models, from mixed to pure states. In 
order to obtain the differential equation corresponding 
to this process, one assumes Markovian behavior, giving 
the evolution of the density operator with reference only 
to its state at present. This Markovian assumption ne- 
glects memory effects: it implies that the time needed for 
the environment to loose the information it received from 
the system is short enough in comparison with the time 
scale of the dynamics we perceive. Then, we are enti- 
tled to regard the information flow in only one direction, 
neglecting any kind of feedback . 

The Lindblad form of the master equation of this 



system-environment model in the Born-Markov approxi- 
mation can be formally written as [T^: 



with formal solution 



p(i) = exp(£t)[p(0)]. 



(1) 



(2) 



where £ stands for the "Lindblandian" operator. In or- 
der to obtain explicit expressions we can trace out the 
environment, which gives 



(3) 



where are the Lindblad operators (/i e [1,...,A^], 
the number A4 depending on the noise model), Hg is the 
system's Hamiltonian and { , } denotes the anticommuta- 
tor. The first two terms of this equation can be regarded 
as the evolution performed by an effective non-Hermitian 
Hamiltonian, Hcs = Hs + iK, with K = —h/2^^L^^L^. 
In fact, we see that 



which reduces to the usual evolution equation for the 
density matrix in the case of i?eff being Hermitian. The 
last term is responsible for the so-called quantum jumps. 
In this context the Lindblad operators are also named 
quantum jump operators. If the initial density matrix is 
in a pure state p(to) = |0(io))(0(io)|i after a time dt 
evolves to the following statistical mixture: 

p{to+dt) = (l-^dp^)|(/)o)(</'o|+^dPMl0M>('^ML (5) 



with the probabilities dp^ defined by 

and the new states by 

, , , [\^iH,adt/h)\<l)[to)) 



1 - dPn 



and 



I0M> = 



\LMh))\ 



(6) 



(7) 



(8) 



Then, the quantum jump picture turns out to be clear; 
with probability dp^ a jump occurs and the system is 
prepared in the state \(f)^). With probability 1 — dp^ 
there are no jumps and the system evolves according to 
the effective Hamiltonian H^s (normalization is included 
also in this case because the evolution is given by a non- 
unitary operator). 
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The numerical method we are going to use in order 
to simulate the master equation is usually known as the 
Monte Carlo Wave function approach • We start from 
a pure state |</)(to)) and at intervals dt, smaller than the 
timescales relevant for the evolution of the density ma- 
trix, we perform the following evaluation. We choose a 
random number e from a uniform distribution in the unit 
interval [0, 1]. If e < dp, where dp — fip^, the system 
jumps to one of the states (to \(f>i) if < e < dpi, 
to 102) if dpi < e < dpi + dp2, and so on). On the other 
hand, if e > dp, the evolution with the non-Hcrmitian 
Hamiltonian 7?eff takes place, ending up in the state |(/)o). 
In both circumstances we renormalize the state. We re- 
peat this process as many times as nsteps = At/dt where 
At is the whole elapsed time during the evolution. Note 
that we must take dt much smaller than the time scales 
relevant for the evolution of the open quantum system 
under investigation. In our simulations, ristops will be 
proportional to the number of quantum gates involved in 
the corresponding protocol. Each realization provides a 
different quantum trajectory and a particular set of them 
(given a choice of the Lindblad operators) is an "unrav- 
eling" of the master equation. It is easy to see that if we 
average over different runs we recover the probabilities 
obtained with the density operator. In fact, given an op- 
erator A, we can write the mean value {A)t ~ Tr[A(o(t)] 
as the average over M trajectories: 

1 ^ 

i—\ 

The advantage of using the quantum trajectories 
method is clear since we need to store a vector of length 
N {N — 2" is the dimension of the Hilbert space, n being 
the number of qubits) rather than a, N x N density ma- 
trix. Moreover, there is also an advantage in computation 
time with respect to density matrix direct calculations. 
We find that a reasonable amount of trajectories (we have 
used 100 < Af < 400 in all calculations, unless otherwise 
mentioned) is needed in order to obtain a satisfactory 
statistical convergence. 

This picture can be formalized by means of the stochas- 
tic Schrodinger equation jl^ 

M0) = -zHMdt - i Y.^LlL^ - {LlL^)m dt 

+ J2 i , -A \^)dN,. (10) 

This is a stochastic nonlinear differential equation, where 
the stochasticity is due to the measurement results: we 
think that the environment is actually measured (as it 
is the case in indirect measurement models) or simpler, 
that the contact of the system with the environment 
produces an effect similar to a continuous measurement 
|l9j| . The nonlinearity is due to the renormalization of 
the state vector. The stochastic differential variables 



dN^ are statistically independent and represent mea- 
surement outcomes. Their ensemble mean is given by 
M[dN^ — {L^^L^) tjjdt. The probability that the vari- 
able dN^ is equal to 1 during a given time step dt is 
{(j)Ljj^Lf^(j)) dt. Therefore, most of the time the variables 
dNf^ are and as a consequence the system evolves con- 
tinuously by means of the non Hermitian effective Hamil- 
tonian. However, when a variable dN^ is equal to 1, the 
corresponding term in equation H10() is the most signif- 
icant. In these cases the quantum jump occurs. Note 
that there are also other possibilities in order to unravel 
the master equation such as the quantum state diffusion 
[Toll20l |. for example. 

B. Quantum operations and quantum trajectories 

There is a close connection between the master equa- 
tion in Lindblad form and the quantum operations theory 
ITsI l and we use this latter formalism as a comparison 
tool for the quantum trajectories calculations. Further- 
more, it will become clear that the evolution of the den- 
sity matrix of the system given by this method can be put 
on the same footing as the stochastic evolution model. In 
fact, the latter constitutes a Monte Carlo simulation of 
the former. 

We write the solution to Eq. (jS)) over an infinitesimal 
time dt as a completely positive map: 

M 

Pit + dt) = $[p{t)] - M^{dt)p{t)Ml{dt), (11) 

where, for /i = 0, we have Mq = 1 — iHcsdt/h and, 
for n > 0, = L^Vdi, satisfying J^^^o H^^t^ = ^ 
to first order in dt. Equation Hll|l is called the Kraus 
representation (or the operator sum representation) of 
the superoperator $ and the operators are known 
as operators elements for the quantum operation $. It 
can be shown that, if the global evolution (system plus 
environment) is unitary, the Kraus operators satisfy the 
completeness relation J2'^o ^t^^-^l^ = 1 Ellll- Note 
that the superoperator $ maps density matrices to den- 
sity matrices, that is p(t+dt) is Hermitian, has unit trace 
and is nonnegative if p{t) satisfies these properties. 

The action of the quantum operation $ can 
be interpreted as p being randomly replaced by 
M^pMl/TT{M^pMl), with probability Tr(M^pMt). 
Equivalently, the set {M^}^=i_... jvi defines a Positive 
Operator Valued Measurement (POVM) with operators 

E^, = MlMf,, that satisfy J2'^o^t^ = ^- The outhned 
process is equivalent to performing a continuous measure- 
ment on the system (or an indirect measurement if the 
environment is actually measured) and shows the close 
connection between the Kraus operators formalism and 
the quantum jumps picture. 

We would like to mention that the quantum opera- 
tions formalism is more general than the master equation 
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approach. The operator sum formulation in differential 
form can be obtained from the master equation, but a 
general quantum process described in terms of an oper- 
ator sum representation needs to be Markovian in order 
to be tractable with a master equation. This opens the 
possibility of studying a wide range of non Markovian 
phenomena using quantum trajectories simulations. 



III. NOISE CHANNELS 

There are several ways to model the interaction of a 
system with the environment. The most common ex- 
amples found in the literature are the amplitude damp- 
ing channel, the phase flip channel, and the depolarizing 
channel "li, Ts"! . In this section, we give a brief descrip- 
tion of two relevant noise models (amplitude damping 
and phase damping) in the single qubit case, and then 
we generalize them to n-qubit systems. 

Dissipation (energy loss) is one of the main features 
present in open quantum systems. Different phenom- 
ena like, for example, spontaneous atomic emission of a 
photon or spin systems approaching equilibrium can be 
modeled by this quantum operation, i.e., the amplitude 
damping. Considering the environment initially in the 
vacuum state |0)e, there is a probability p that the ex- 
cited state of the system decays and that the state of the 
environment changes from the vacuum to 

|0),|0)e ^ |0),|0)e, 

|l),|0)e ^ v/r^|l).,|0)e + VP|0)s|l)e, (12) 

where |0)s stands for the spin down (ground) state and 
\l)s for the spin up (excited) state of the system, and 
the indexes s and e denote the quantum states of the 
system and of the environment, respectively. Tracing out 
the environment the corresponding Kraus operators are 
obtained: 



with Kraus operators 



Mo 







1 







Ml 



VP 




(13) 



The operator Mi is responsible for the quantum jumps 
and Mq for the continuous evolution. For the case of in- 
finitesimal evolution operators (see Sec. a repeated 
application of this noise channel gives an exp onential de- 
cay law of the population of the state |1) [l3- Therefore, 
this evolution drives any intial (pure or mixed) state of 
the qubit to the pure state |0). Note that here and in the 
following the matrix representations of the single-qubit 
Kraus operators are written in the {|0), |1)} basis. 

We also consider the phase flip channel (which is equiv- 
alent to the phase damping channel jl^) given by the 
following model 



|0)s|0), 
|1>.|0). 



yi^|0)3|0),+^|0),|l)e, 
v/T^|l),|0)e-VP|l>.|l>e, (14) 



(15) 



This kind of noise can be thought as describing quantum 
information loss, in contrast with the previous model, 
which describes energy loss. This is a purely quantum 
mechanical process and has been extensively studied in 
the context of quantum to classical correspondence ITqI . 

In order to generalize the single qubit amplitude damp- 
ing process to many qubits we will follow two different 
points of view. In the first case we assume that a single 
damping probability describes the action of the environ- 
ment, irrespective of the internal many-body state of the 
system. In the second approach, we assume that each 
qubit has its own interaction with the environment, in- 
dependently of the other qubits. This makes the damping 
probability grow with the number of qubits that can per- 
form the transition |1) |0). Both models assume that 
only one qubit of the system can decay at a time. 

In the first case we have a probability p for the 
system to perform one of the possible transitions, each of 
them being equally likely. This can be illustrated with a 
two- qubit example: 



|00),|00), 
|01).|00), 
|10).|00), 
|11).|00), 



|00),|00)e, 

v/r^ioi),|oo)e 



yT^|io).|oo), 
v/T^|ii).|oo), 

(|10),|01)e + |01),|10)e). 



Vp|oo).|oi), 
Vp|oo).|io), 

Vp72 



(16) 



A n-qubit state |i„_i...«o) (ii = 0,1, with < 
j < n — 1) decays in the interval dt with a probabil- 
ity p — Tdt/h. After this infinitesimal time, the possi- 
ble states of the system are those in which the damping 
|1) |0) has occurred in one of the qubits, the damping 
probability being the same for all the qubits. We provide 
a compact expression for the Kraus operator s M,, in our 
model, following the formulation given in SecEHfor the 
infinitesimal evolution. The matrix elements for the fcth 
operator in the computational basis |i) = |i„_ 
with i = X^zLo are given by 



■to) 



, forj>2(''-i). 



- n - 2(^-1) , N 







otherwise. 



where the superscript ^"^ in M^" underlines the fact that 
we are dealing with n-qubit Kraus operators. There are n 
operators M^"^ (m = Ij ■ • • i 't^); where the index singles 
out which qubit undergoes the transition |1) — > |0). The 
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A/q"'' operator is 



(»)l _ 



1 - 







for i = j = 0, 
^ for i = j 7^ 0, 

otherwise. 



(18) 



Note that, if we replace in the above definition the square 
root by its first order approximation, we arrive at the 
same expression for the action of the effective Hamilto- 
nian given in Sec. (see the first term in the right hand 
side of Eq. © and Eq. O). 

For example, starting from the four-qubit pure state 
p{to) = |1011)(1011|, the action of the generalized am- 
plitude damping channel leads, after a time dt, to the 
statistical mixture 



p{to + dt) = 



1-I|!')|1011)(1011| 



^(|0011)(0011| 



|1001)(1001| + |1010)(1010|). 



(19) 



This situation can be described as branching process or 
a cascade in the population of different classes of states 
(see, e.g., Ref. |21|). A class in this model is naturally 
defined as the collection of all the states of the system 
having the same number of qubits in the "up" state (i.e., 
of |1) states). The first class is the one initially populated 
with probability Wq. The second class, with associated 
probability Wi, corresponds to the states that are ob- 
tained from the initial state after the damping of one 
qubit. The third class, with probability W2, includes the 
states that can be reached by the damping of one qubit 
in any of the states of the second class, and so on until 
reaching the "ground state" |0 . . . 0) of the system. This 
last state is also the only state that resides in the last 
class. This class is populated with probability Wm, the 
number m being the maximum number of up spins in 
the initial state (if we start from a state \in-i . . . *o) of 
the computational basis, then m — Yl"=o Note that 
m < n. The set of probabilities Wm can be obtained by 
means of the following differential equations: 

dWo r 



dWk 
dt 



dW„ 
dt 



{Wk-i - Wk), 



r 

n 



(20) 



with solutions 



w„ 



m— 1 
fc=0 



(21) 




FIG. 1: Probabilities Wfc for each class of states in terms of the 
dimensionless damping rate 7 = Tdt/h, in a system with n = 6 
qubits, obtained after evolving the initial state |1 . . . 1) up to time 
2ndt, under the noise model I17i . Solid lines correspond to the 
exact formulae 1211 . and the numbers close to each curve indicate 
the corresponding class. Diamonds stand for quantum trajectories 
simulations (error bars are not shown since they are smaller than 
the size of the symbols). 



In Fig. n we show a comparison among the probabilities 
for each class obtained with quantum trajectories simu- 
lations and the theoretical formulae, in the case of n = 6. 
We have evolved the initial state |2" — 1) = |1 . . . 1) up 
to time 2n dt, for different values of the dimensionless 
damping rate 7 = Tdt/h. This evolution is purely dis- 
sipative, without considering any other kind of system's 
dynamics. As can be seen from Fig. ^ we have a very 
good agreement between quantum trajectories numerical 
simulations and the theoretical predictions of the model. 

The other generalization for the amplitude damping 
differs from the previous one in that the decay probabil- 
ity is now the same for each single qubit process, inde- 
pendently of the state of the system. Then, the decay 
probability for a states of the computational basis is pro- 
portional to the number of qubits in the up state. This 
can be illustrated in the two-qubit case. In this example, 
the noise channel is described by the following unitary 
evolution formulae: 

|00),|00)e ^ |00),|00)e, 

|01),|00)e ^ yT^|01),|00)e + ^|00),|01)e, 

|10),|00)e ^ v/^|10).|00>e + \^|00),|10)e, 

|ll),|00)e ^l-2p\n)sme + VP 

(|10),|01)e + |01),|10)e). (22) 

(n) 

The Kraus operators M/i for this model are given by 
the n-factor tensor product 

a (g) ...» Ml «... (g) 1, (23) 
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where Mi is given by Eq. H13() and ^ {ji = 1, . . . , n) co- 
incides with the position of Mi in that is /i singles 
out the qubit which decays from |1) to |0). In the compu- 
tational basis, the matrix representation of the operator 
Mq""* is given by 



(«)i 



1 ioY i — j ~ 0, 

l-Er=o'*;^ forz = j^O, 



(24) 



otherwise. 



The evolution of the pure state p(to) — |1011)(1011| is 
different from what obtained in the previous many-qubit 
damping model (see Eq. (|19|l . We have 



-dt) = ( 1 - |1011>(1011|+^|^(|0011)(0011| 



|1001)(1001| + |1010)(1010|). 



(25) 



The cascade in the population of the different classes 
Wk is ruled by the following set of differential equations: 

rfWo r 



dWk 
dt 



dW„ 
dt 



{nk-i Wk-i - rik Wk), 



n 



-1 W^m-l = -^Wm-l, 

n 



(26) 



where is the number of qubits in the "up" state for the 
class k. Taking into account that Uk = ng — k = m — k, 
the solutions are 



n\ ^ (-l)(fe-') r n,,rt 



nfc! ^ i! (fc-i)! 

m— 1 



(27) 



k=0 



In Fig. 12 we take the initial state |2" — 1) evolving it 
for a time n dt, for different values of the dimensionless 
damping rate 7. Again, only dissipation is considered and 
we have found a very good agreement with our analytical 
predictions. 

With this same idea we generalize the phase flip chan- 
nel. The Kraus operators are given by the n-factor tensor 
product of Eq. (|23|l . where now ^ (/i = 1, . . . , n) coincides 
with the position of the matrix Mi of Eq. IjlSI) , instead 
of Eq. . The matrix representation of the Kraus op- 
erator Mg"'' in the computational basis is 



l-n^ fori = j. 



otherwise. 



(28) 




FIG. 2: Same as in Fig. but for the noise model 1231 and after 
evolution up to time ndt. 



We study the stability of a given state vector subjected 
to this decoherence channel. Namely, we compute the 
fidelity F = Tr[/5o/o(i)] of a n-qubit system, for a ran- 
dom initial state {ipo), as a function of 7 = Tdt/h 
{Pq — |V'o)(^o| is the density matrix for the initial state 
and p{t) is the density matrix of the system at time t). 
All components of the initial random state vector lipo) 
are taken to be random complex numbers of modulus 
l/\/2". A combinatorial calculation provides an exact 
closed expression for the fidelity in this case. We obtain 



1 I " 



1 



{n-i)\ 



exp 



2i rt 



(29) 



Note that this formula does not give a simple exponential 
fidelity decay, but a superposition of exponential decays 
with different rates, since the decay rate depends on the 
number of up spins in the different states of the compu- 
tational basis. 

The agreement of the quantum trajectories simulations 
with this theoretical formula is shown in Fig.|3| The semi 
log version of the theoretical curve is shown in the inset, 
where we plot F — Foo , Foe — 1 /2" being the asymptotic 
value of fidelity at t = 00 or at F = 00. We note that the 
asymptotic decay of F = F — Foo takes place with the 
lowest decay rate Fi = 2F/?i. of Eq. it^ . 

We note that analytical formulae for fidelity decay can 
also be found for other special initial conditions. For 
instance, Eq. (|29|l remains valid when the initial state 
is a random superposition of computational basis states 
with up to ni spins up, provided that we replace n with 
ni in this equation. 

In the following, we will apply the amplitude damping 
models to study the fidelity of a generalized teleporta- 
tion protocol and the phase flip channel for the case of a 
quantum computer implementation of the baker's map. 
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0.04 



FIG. 3: Fidelity of a random initial state in terms of the dimension- 
less decay rate 7 = Vdt/h, in a system with n = 6 qubits, subjected 
to the generalized phase flip channel up to time 2n^dt. The very 
good agreement between quantum trajectories simulations (circles 
with error bars) and theoretical predictions (solid line) is clearly 
seen. Here and in the following figures the error bars give the size 
of the statistical error. The inset shows the semi log version of the 
theoretical curve plot for the modified fidelity F = F — Faa ■ 



Bob 



System 



Alice 



87 654321 
Environment 

FIG. 4: Schematic drawing of the teleportation procedure that 
we study in the text. Alice sends one of the qubits of her pure 
Bell state to Bob. Meanwhile there is dissipation induced by the 
environment. In the figure, there is a chain of n = 9 qubits and the 
third of the ra — 2 swap gates required by this quantum protocol 
has been applied. Qubit |»/)) has to be teleported. 



IV. QUANTUM TELEPORTATION 

Recently, there have been several publications focused 
on the investigation of fidelity of teleportation in the pres- 
ence of a noisy environment Q2^ 23., .24, 2M- We study 
this problem in the situation presented in jl^l, where a 
model of quantum teleportation through a noisy chain 
of qubits has been used. A schematic drawing of this 
quantum protocol is shown in Fig. 0] 

Let us first recall this protocol in the ideal case, with- 
out environmental effects. We consider a chain of n 
qubits, and assume that Alice can access the qubits lo- 
cated at one end of the chain. Bob those at the other 



end. Initially Alice owns an EPR pair (for instance we 
take the Bell state |(/)+) = (|00) + \ll))/y/2), while the 
remaining n — 2 qubits are in a pure state. Thus, the 
global initial state of the chain is given by 



E 



.,Z2 l^n — 1 • 



■*2)®^(|00) 



111)), (30) 



where ik = 0,1 denotes the up or down state of the qubit 
k. In order to deliver one of the qubits of the EPR pair 
to Bob, we implement a protocol consisting of n — 2 swap 
gates that exchange the states of pairs of qubits: 



E 



E 



V2 



V2 



. ^200) 



.0^20) + \i 



ri-l • • 



.^211)) 



H2I)) 



E 



V2 



-(|0«r, 



..j20) + |li„-l...i2l)). 

(31) 

After that, Alice and Bob share an EPR pair, and there- 
fore an unknown state of a qubit (j?/') = a|0) + h\l)) can 
be transferred from Alice to Bob by means of the stan- 
dard teleportation protocol .13]. In this work, we take 
random coefficients Ci^_-^^^...^i2, that is they have ampli- 
tudes of the order of l/-\/2"~2 (to assure wave function 
normalization) and random phases. This ergodic hypoth- 
esis models the transmission of a qubit through a chaotic 
quantum channel. 

We assume that our quantum protocol is implemented 
by a sequence of instantaneous and perfect swap gates, 
separated by a time interval t, during which the corre- 
sponding noise channels introduce errors. This means 
that, using quantum trajectories, dissipation is imple- 
mented by means of "infinitesimal" Kraus operators (see 
Sec. IIII|I . following the quantum jump numerical proce- 
dure outlined in Sec. Ill AI We also assume that the only 
effect of the system's Hamiltonian Hs is to generate the 
swap gates. 

in) 

Let us call pj. the density matrix of the whole chain 
of n qubits after k swap gates. Since the evolution of the 
density matrix in a time step dt is given, in the Kraus 
representation, by Eq. Hll|) . we can write the evolution 

from p^^li to /9^"^ as follows: 



„(") _ TTk,k+l 



M 



J2 M^iHdt) ■ ■ ■ M^-Hdt) 

./ii,...,/ifc=0 

fe,fc+it 



(32) 



where J7g4 is the swap operator that exchanges the states 
of the qubits i and j and k = r/dt is the number of 
quantum noise operations between two consecutive swap 
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gates. After n — 2 swap gates, we obtain the final state 

(") (") 

Pf = Pn-2- 

After computing the evolution of the initial state of the 
chain of n qubits up to time At = (n — 2)t , the standard 

m 



teleportation protocol is implemented 
of teleportation is defined by 



F 



\Pb\ 



The fidelity 



(33) 



where — a\Q) +b\l) is the state to be teleported, and 
Pb is the density matrix of Bob's qubit at the end of the 
teleportation protocol, obtained after tracing over all the 

other qubits of the chain: pB = Tro,...,n-2[pj"']- 

In the quantum trajectories method, we compute the 
fidelity as 



i—l 



(34) 



where {pB)i is the reduced density matrix of Bob's qubit, 
obtained from the wave vector of the trajectory i at the 
end of the quantum protocol. If the final state of the 

chain is |0) = X]j=o^ Q^ili) (^^ arbitrary state), the state 
of the whole system is 

+ {aa2f - ba2j'+i)\i')\(j)^). 
+ {ba2f + aQ;2j'+i)|/)|V'^) 
+ (6q2j' - aQ;2j'+i)|j')|V'")], (35) 

In the previous expression we have used the four (max- 
imally entangled) Bell states (|0±) = (|00) ± |ll))/\/2 
and = (|01) ± |10))-\/2), corresponding to the two 

least significant bits (the first bit of the chain and |-0)). 
Then, as required by the teleportation protocol, we per- 
form a Bell measurement whose result determines one 
out of possible unitary transformations acting on Bob's 
qubit. Owing to quantum noise, Bob's qubit is entangled 
with the rest of the chain. Therefore, we must trace over 
these qubits to obtain the reduced density matrix ps de- 
scribing the state of Bob's qubit. If we measure \(t>'^) for 
instance, and define oiji ~ l/\/2 (a a2ji + b a2j'+i), we 
arrive at the following expression: 




where D 



2"-2 and C = 



2D-1 
J'=0 



is a normal- 



ization constant. Then, the fidelity becomes 



F= C{\a\^Y.T=^l\ar? + \b?Y:7jj^ 



+ E,^Io2Re(a,v5;,+^a*&)). (37) 
In the special case a = b = 1/ V2, it reduces to 



1 



D-l 

i^=- + C^Re(5ya*,+^) 
j'=o 



(38) 



Similar expressions are obtained when the Bell measure- 
ment gives outcomes or 

Let us now briefly discuss the case in which we directly 
evolve the density matrix of the whole system by means of 
the master equation Q. After the sequence of swap gates 
and quantum noise operations, we obtain the final density 

(n) 

matrix pj for the rt-qubit chain. The Bell measurement 
is performed by the operator 



Tr[M^^{p\"'> (g>p^)Ml 



(39) 



where p^ = |'0)('0| is the density matrix describing 
the state of the qubit to be teleported and = 
a^""^^ (g) 10+) (0+1 (a^""^) is the identity for the remain- 
ing n — 1 qubits of the chain). Then, we take the partial 
trace and obtain the fidelity F — Tr(pBPV') ~ (^IPbIV')- 
It is straightforward to see that the above procedure is 
equivalent to take the partial trace (over the mediating 
qubits of the chain) first and then proceed with low di- 
mensionality calculations. 

We have investigated the effect of the two amplitude 
damping models described in Sec. IIIII on the fidelity of 
quantum teleportation. In Fig. we show the results of 
our numerical simulations for the special case in which 
the state to be teleported is lip) — (|0) + |l))/-\/2. In this 
case, in the limit of infinite chain (n —^ oo) or of large 
damping rate, the density matrix pe describing the state 
of Bob's qubit becomes pB = |0)(0|. Thus, the asymp- 
totic value of fidelity is given by Fao = 1/2 and we plot 
the values of F = -F — Foo, for the noise models H17|l 
and (|23|l . corresponding to the inset and the main figure 
in Fig. [SJ respectively. For both cases we have checked 
the accuracy of the quantum trajectory simulations by 
reproducing the results with direct density matrix calcu- 
lations. This was possible only up to n = 10 qubits. The 
exponential decay of fidelity in the case in which quan- 
tum noise is modeled by Eq. H23|) is in agreement with 
the theoretical formula 



11 /Ft 
F = - + - exp 

2 2 ^ I n 



i + iexp(-7A:), (40) 



where 7 = Tt/Ti is the dimensionless damping rate and 
k — t/r = n — 2 measures the time in units of quantum 
(swap) gates. To derive this theoretical formula, we ob- 
serve that this quantum noise model does not generate 
entanglement between the two qubits of the Bell pair and 
the others qubit of the chain. Therefore, it is sufficient 
to study the evolution of the Bell state |0+)(0+| under 
the noise model (|23|l . This evolution takes place inside 
a two-qubit Hilbert space of dimension 4 and its exact 
solution is given by Eq. 100)). On the contrary, the quan- 
tum noise model ((T7I) entangles these two qubits with 
the rest of the chain. In this case, the fidelity decay is 
not exponential. Unfortunately, we could not provide an 
analytical derivation of F{t) in this case. 



9 



10" 



10" 



10" 



10" 



10" 



10" 




0.4 0.5 



FIG. 5: Fidelity F = F-F^ (Foo = 1/2) of the teleportation 
of the state \ip) = (|0) + \1))/V2 as a function of the dimen- 
sionless damping rate 7 — Tr/Ti, for the amphtude damping 
model 1231 . Circles and squares are the results of the quantum 
trajectories calculations for chains with n = 10 and n = 20 
qubits, respectively. Triangles give the results of the density 
matrix calculations at n = 10. Straight lines correspond to 
the theoretical result of Eq. I|40^ . Inset: the same but for the 
noise model 1171 . 



V. THE QUANTUM BAKER'S MAP 

The quantum algorithm for the simulation of the quan- 
tum baker's map has been proposed by Schack and 
recently implemented by means of a three-qubit NMR- 
based quantum processor where the fidelity of the 
quantum computation of the baker's map in one map 
step was measured. 

The baker's transformation is one of the prototype 
models of classical and quantum chaos [l^. It maps the 
unit square < g,p < 1 onto itself according to 



qk+i 
Pk+i 



2qk - [2qk], 
(pk + [2qfe])/2, 



(41) 



where [x] stands for the integer part of x and the index k 
denotes the number of map iterations. The action of this 
map corresponds to compressing the unit square in the p 
direction and stretching it in the q direction, then cutting 
it along the p direction, and finally stacking one piece on 
top of the other (similarly to the way a baker kneads 
dough). Note that map H41|l is area preserving. The 
baker's map is a paradigmatic model of classical chaos. 
Indeed, it exhibits sensitive dependence on initial condi- 
tions, which is the distinctive feature of classical chaos: 
any small error in determining the initial conditions is 
amplified exponentially in time. In other words, two 
nearby trajectories separate exponentially, with a rate 
given by the maximum Lyapunov exponent A = In 2. 
The baker's map can be quantized following |2^. We 



introduce the position {q) and momentum {p) opera- 
tors, and denote the eigenstates of these operators by 
\qj) and \pk), respectively. The corresponding eigen- 
values are given by qj = j/N and pk — k/N, with 
j, k — 0, . . . , N — 1, N being the dimension of the Hilbert 
space. Note that, to fit N levels onto the unit square, we 
must set 2Trh = 1/N. Therefore, the effective Planck's 
constant of the system is hcs cx 1/A^. We take N = 2", 
where n is the number of qubits used to simulate the 
quantum baker's map on a quantum computer. Note 
that ?icff drops exponentially with the number of qubits 
and therefore the semiclassical regime heg <C 1 can be 
reached with a small number of qubits. The transfor- 
mation between the position basis {\qo), . . . , \qN-i)} and 
the momentum basis {|po); ■ ■ ■ 1 \pn-i)} is performed by 
means of a discrete Fourier transform defined by the 
matrix elements 



{qk\Fn\qj) 



1 



/ 2'Kikj\ 



(42) 



It can be shown 26] that the quantized baker's map may 
be defined by the transformation 



IV', 



fc+1/ 











Hk), (43) 



where \^pk) denotes the wave vector of the system after k 
map steps, the matrix elements are to be understood rel- 
ative to the position basis {l^j}} and -Fn-i is the discrete 
Fourier transform, defined by Eq. (|42|l . 

Since the discrete Fourier transform can be calculated 
on a quantum computer using 0{n^) elementary gates 



(see, e.g., Ref. the simulation of one ste- 

baker's map requires 0(n^ = (logiV)^) gates 



of the 
(more 



precisely, we need (n — 1) controlled phase-shift gates 
and 2n — 1 Hadamard gates). Therefore, it is exponen- 
tially faster than the best known classical computation, 
which is based on the fast Fourier transform and requires 
O(A^logiV) gates. 

In this section, we investigate the fidelity of the quan- 
tum computation of the baker's map in the presence of 
quantum noise. We consider the phase flip noise channel, 
generalized to the n-qubit case as discussed in Sec. IIIII 
We take an initial state jV'o) with amplitudes of the order 
of l/-\/2" and random phases. We perform the forward 
evolution of the baker's map up to time k (this evolution 
is driven, in the noiseless case, by the unitary operator 
B'', with B given in Ea. l43|l . followed by the /c-stcp back- 
ward evolution (represented by the operator (B^)''). Due 
to quantum noise, the initial state |V'o) is not exactly re- 
covered and the final state of the system is described by a 
density matrix pf. The fidelity of quantum computation 
is given by F = (^/iqI/O/IV'o)- 

We are able to work out the following theoretical for- 
mula for the decay of fidelity induced by phase flip noise: 



F — exp {—n^Ng) — exp{—2jn^k), 



(44) 



where 7 = Tr/h is the dimensionless decay rate (again, 
T denotes the time interval between elementary quantum 
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gates) and Ng = 2-n?k is the total number of elemen- 
tary quantum gates required to implement the k steps 
forward evolution of the baker's map, followed by k step 
backward. To derive this formula, we first of all note that 
Eq. H29I) . obtained in the absence of any quantum gate 
operation, gives, at short times, the exponential decay 



10" 



<V>f 



(45) 



Here the decay rate < F > is obtained after averaging all 
the decay rates that appear in Eq. (|^ : 



2" ^ 



i=0 







2iT = nr. 



(46) 



The effect of the chaotic dynamics of the baker's map is 
to induce a fast decay of correlations, so that any mem- 
ory of the initial state is rapidly lost and the fidelity 
decay remains exponential also at long times. Therefore, 
the condition for the validity of formula (|44|l is that the 
randomization process introduced by chaotic dynamics 
takes place in a time scale shorter than the time scale for 
fidelity decay. 

We can determine from Eq. (|44|l the time scale up to 
which a reliable quantum computation of the baker's map 
evolution in the presence of the phase flip noise channel 
is possible even without quantum error correction. The 
time scale kf at which F drops below some constant A 
(for instance, A = 0.9) is given by 



kf = - 



liiA 



(47) 



The total number of gates that can be implemented up 
to this time scale is given by 



7n 



(48) 



Our theoretical expectations are confirmed by the nu- 
merical data of Figs. and [7| In Fig. El we show the 
fidelity after one map step, as a function of the dimen- 
sionless damping rate 7 — Tt/Ti. The numerical data 
of this figure show that the fidelity drops exponentially 
with 7, in excellent agreement with Eq. H44|l . We point 
out that our theory predicts not only the exponential fi- 
delity decay but also the numerical value of the decay 
rate. Finally, we show in Fig. Qthe number kf oi for- 
ward/backward map steps required to reach a fixed value 
A of the fidelity {F — A — 0.9), as a function of the num- 
ber n of qubits. This time scale decays as a power law, 
kf = C/n^, in agreement with Eq. (|44|l . Again, our the- 
ory predicts also the value of the prefactor C, in excellent 
agreement with numerical data. 

As we have seen above, the number of gates that can 
be reliably implemented without quantum error correc- 
tion drops only polynomially with the number of qubits, 
(iVg)/ cx 1/n (see Eq. We would Hke to stress that 
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FIG. 6: Semi log plot of the fidelity as a function of the 
dimensionless decay rate 7, for the baker's map after one map 
step in the presence of the phase flip channel described in the 
text. Circles and squares correspond to quantum trajectories 
simulations with A/" = 500 trajectories, at n = 10 and n = 20 
qubits, respectively. Triangles give the results obtained by 
direct computation of the density matrix evolution at n = 10. 
Solid lines stand for the theoretical prediction of Eq. 11441 . 
namely ^(7) = exp(— 2771^). 
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FIG. 7: Logarithmic plot of the time scale kf (measured in 
number of steps of the baker's map) required to reach a fidelity 
value of F — 0.9, as a function of the number n of qubits, for 
7 ~ 6.47 X 10~^ (this value has been chosen by demanding 
that F — 0.9 after one forward/backward step of the baker's 
map ai n — 21). The straight line gives the theoretical curve 
T{n) = Cln\ with C = - ln(0.9)/(27) « 8.14. 
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this dependence should remain vahd also in other envi- 
ronment models that allow only one qubit at a time to 
perform a transition, like the other noise channels des- 
cussed in Sec. IIIII Furthermore, we note that the time 
scales for fidelity decay derived in this section and con- 
firmed in the baker's map model, are expected to be valid 
for any quantum algorithm simulating dynamical systems 
in the regime of quantum chaos. 

VI. CONCLUSION 

We have studied two quantum protocols, i.e. a tele- 
portation scheme through a chain of qubits and a quan- 
tum algorithm for the quantum baker's map. We have 
modeled different sorts of environments in order to get a 
deeper insight of the stability of quantum computation 
when different interactions with the environment taken 
into account. Two kinds of generalized amplitude damp- 
ing models have been considered for the teleportation 



scheme, giving very different behaviors. After a theoreti- 
cal analysis we were able to understand the origin of these 
differences. This reveals the importance of the details 
of the environmental models in assessing the operability 
bounds of quantum processors. In the case of the baker's 
map simulation, wc have chosen the phase flip type of 
noise and we could verify our theoretical predictions for 
fidelity decay. The results of this paper show that quan- 
tum trajectories are a very valuable tool for simulating 
noise processes in quantum information protocols with a 
high degree of efficiency. 
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